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ABSTRACT 

We present here dust models for the Small Magellanic Cloud bar calculated for the first time with the 
regularization approach. A simple mixture of the core-mantle and/or composite grains (mostly made 
from silicates and organic refractory) together with silicon nanoparticles is consistent with the following: 
1) the observed extinction toward AzV 398 (a typical SMC bar sightline), 2) the elemental abundances, 
and 3) the strength of the interstellar radiation field. We predict the expected albedo and asymmetry 
parameter of the models, which can be tested in future observations. The proposed dust models can also 
be tested by looking for the expected extended red emission. 

Subject headings: dust, extinction - galaxies: individual (Small Magellanic Cloud) - galaxies: ISM - 
Magellanic Clouds 



1. INTRODUCTION 

The Small Magellanic Cloud (SMC), one of our nearest 
neighbours, contains considerably less heavy elements and 
dust than the Galaxy (Bouchet et al. 1985; Welty et al. 
1997). In addition, dust in the SMC see ms to be di fferen t 

' Pei pgg ). 
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rom dust in the Galaxy (Prevot et al. 
For example, a typical extinction curve of the SMC is al- 
most linear with inverse wavelength and does not s how a 
presence of the UV bu mp at 2175 A (Prevot et al. [1984| ; 
Thompson et al. 1988). Recently, Rodrigues et al. ( |1997| ) 
have measured the linear polarization for a sample of the 
SMC stars in the optical and found that the wavelength 
of maximum polarization is generally smaller than that in 
the Galaxy. The low metallicity of the SMC implies that it 
should be at an early stage of its chemical evolution, thus 
resembling in this respect galaxies at high redshifts. A 
strong support to this view has come from the discovery 
that the dust in starburst galaxies, apparently the only 
type of galaxies found so far for redshifts z > 2.5, has 
an extinction law remarka bely si milar to that in the SMC 
(Gordon, Calzetti, & Witt 09T\). 



Several attempts have been made to model the d ust in 
the SMC. Bromage & Nandy ( |l983| ) and Pei ( |l992| ) mod- 
elled some SMC mean extinction curves using the dust 
mixture of spherical graphite and silicate grains with a 
powe r-law size distribution (Mathis, Rump l & Nordsieck 
1977| , h ereaft er MRN; Drain e & Lee |1984[ ). Bromage & 
Nandy ( |1983| ) and Pei ( |1992| ) have shown that an MRN- 
like mixture with a lower fractional mass of graphite grains 
in comparison to silicate grains, and with other parameters 
as in the Gal axy, c an satisfactorily explain the SMC ex- 
tinction. Pei (1992) even succeeded to fit the SMC extinc- 
tion law with silicate grains alone. Rodrigues et al. (1997) 



made model fits to both extinction and polarization for two 
stars in the SMC, AzV 398 and AzV 456, representing lines 
of sight with different properties. The authors have found 
that the mixture of bare silicate and amorphous carbon, 
or graphite spheres together with silicate cylinders with an 
MRN-like power-law size distribution explains quite well 
both the extinction and the polarization. 



However, the SMC dust models proposed so far are 
too simplistic concerning both the choice of grain con- 
stituents and grain-size distri butions. Recently, Mathis 
(p9fj), Li k G reenberg ( pfjj ) and Zubko, Krelowski, & 
Wegner (1998) have presented more sophisticated models 
of Galactic dust using core-mantle, multilayer and com- 
posite grains which are much more physically justified then 
the models previ ously considered. Note that the models 
by Zubko et al. 



w ere ca lculated with the regular- 
ization approach (Zubko 1997) which is a very efficient 



method capable of deriving optimum and unique size dis- 
tributions in a general form for any predefined mixture of 
grains by simultaneously fitting the extinction curve, the 
elemental abundances and the mass fraction constraints. 
The uniqueness of the grain size distributions follows from 
the mathematical nature of the problem when we need to 
solve a Fredholm integral equation of the first order, being 
a typical ill-posed problem. The regularization approach 
reduces this problem to minimization of a strongly con- 
vex quadratic functional. The latter problem was strictly 
proved to have a single solution. S ee for more detai ls, e.g . 
Groetsch ( pgg ), Tikhonov et al. ( p9Cj ) or Zubko ( |l997| ). 
This method can be expanded to allow one to deduce the 
uncertainty in the solution based on the data uncertainty, 
but this is beyond the scope of this letter (Zubko 1999, in 
preparation) . 



Recently, Gordon & Clayton (1998) have derived the ex- 



tinction curves for four SMC stars with several improve- 
ments in comparison to previous studies: higher S/N IUE 
spectra, and a more careful choice of the pairs of red- 
dened and comparison stars. The sightlincs toward three 
stars: AzV 18, 214, and 398, located in the SMC bar pass 
through the regions of active star formation and exhibit 
similar extinction law. The sightline toward the star AzV 
456, located in the SMC wing, passes through a much more 
quiescent region of star formation and shows a Galaxy- like 
extinction with the 2175 A UV bump. The purpose of the 
present Letter is to report the first SMC dust models cal- 
culated with the regularization approach to the new high 
quality extinction curves. We modelled the extinction to- 
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ward the star AzV 398 which is t hough t to be a typical 



SM C bar sightline (Rodrigues et al. 1997; Gordon & Clay- 



ton 



199S). The results of the study which includes all the 



four stars will be presented in a forthcoming paper. 



2. EMPIRICAL DATA 
We transform the extinction curve for AzV 398, de- 



rived by Gordon & Clayton (1998) in the standard form: 
E(X) = [t(X)-t(V)]/[t(B)-^ 
section per H atom: 



to the extinction cross 



Ap 



Ap 



Rv] 



(1) 



where r is the optical thickness, E(B — V) is the 
B — Ucolour excess, Ry is the total-to-selective extinction 
ratio, and Ah is the column number density of hydrogen. 
For Ah w e take the value 1.5xl0 22 cm -2 from Bouchet et 
al. ( 1985Q , corresponding to atomic hydrogen. Since the 
sightline toward Az V 398 is associated with an H II region 
(Gordon & Clayton 1998] ), it is likely that the contribution 
of molecular hydrogen to Ah is neglegibly small (for sight- 
lines not passing t hroug h the SMC bar this may not be the 
case, see Lequeux 1994 ). The val ue for E(B — V)— 0.37 was 
also taken from Bo uchet et al. ( 1985| ) and Rv=2.87 from 
Gordon & Clayton fll99S| ). 

The elemental abundances (gas + dust), cu rrentl y 
adopted for the SMC, were taken from Welty et al. ( |1997| ): 
C/H=46 p.p.m. (atoms per 10 6 H atoms) or 7.66±0.13 
dex, O/H=107 p.p.m. or 8.03±0.10 dex, Si/H=10 p.p.m. 
or 7.00±0.18 dex, Mg/H=9.1 p.p.m. or 6.96±0.12 dex, and 
Fe/H=6.6 p.p.m. or 6.82±0.13 dex. Note that these values 
are 2-5 times less than the respective Gal actic abundances 
foll owing recent revision (Snow & Witt 1996; Cardelli et 



al. 1996). Since we have no information on the amounts of 



elements in dust and gas separately in the SMC, we simply 
assume in this study that as in the Galaxy 42% of carbon 
(~20 p.p.m.), 37% of oxygen (~40 p.p.m.), and all silicon, 
magn esium and iro n are locked up in dust (Cardelli et al. 



1996; Zubko et al. 1998] ). The actual amount of elements 



locked up in dust is uncertain, but one may expect even 
lower amounts of elements in dust because the SMC is less 

chemically processed than our Galaxy. 

Rec ently, Witt, Gordon, & Furton ( 1998] ) and Ledoux 
et al. (1998) proposed that the silicon nanoparticles might 
be the source of the extended re d emi ssion (ERE) in our 
Galaxy. Zubko, Smith, & Witt ( |1999| ) have modeled the 
mean Galactic extinction curve with the silicon nanopar- 
ticles involved and have shown that this hypothesis is con- 
sistent with the available data on extinction, elemental 
abundances an d the ERE. On the other hand, Perrin , Dar- 
bon, & Sivan fll995| ) and Darbon, Perrin, & Sivan fll998| ) 
have revealed ERE in extragalactic objects showing active 
star formation: the starburst galaxy M82 and the nebula 
30 Doradus in the Large Magellanic Cloud, respectively. 
Since the sightline toward AzV 398 passes through a star- 
forming region, we may expect to observe ERE from there 
as well. We thus included the sili con n anoparticles in our 
modeling. As in Zubko et al. (199E), we used the sili- 
con core-Si02 mantle model of silicon nanoparticles with 
optica l constants of nanosized silicon from Koshida et al. 
( 1993| ). We also included in present study the grain con- 



stituents (graphite, silicate, SiC, organic refractory, amor- 
phous carbon, water ice and others) and res pectiv e optical 
constants previously used by Zubko et al. ( 1998| ). 



3. MODELS OF EXTINCTION 

We performed extensive work on modeling the extinc- 
tion curve for AzV 398, searching for the physically rea- 
sonable mixtures of dust constituents. Our goal was to 
find the models which would simultaneously fit the extinc- 
tion, consume the allowed amounts of chemical elements 
and include silicon nanopa rticle ss (by analogy with the 



Galactic case, Zubko et al. 1999). We report in this Let- 



ter three simple models which fulfil all the above require- 
ments. The results are presented in Figs [l|-|2| and Table |l| 
The model grains are mostly made up of two species: sil- 
icate (MgFeSi04) and organic refractory residue, which 
coexist either in core(silicate)-mantle(organic refractory) 
or in spherical porous composite grains with the latter 
containing also small amounts of amorphous carbon. The 
total mass fraction of silicate + organic refractory is about 
0.9. The other important model component is the silicon 
nanoparticles, which are found to have a mass fraction of 
about 0.07-0.085. 

The fact that organic refractory is among the major 
grain constituents is in good agreement with the expec- 
tations that the interstellar radiation field (ISRF) in the 
SMC is stronger than in the Galaxy (Lequeux et al. [1994 ) 
since the icy mantles on silicate grains formed in molecular 
clouds can be processed by t he UV radiation into organic 
refractory (Greenberg & Li 1996). Note especially that 
our attempts to include silicate and SiC grains coated by 
either amorphous carbon or water ice mantles and also 
bare carbonaceous (graphite, amorphous carbon), silicate 
and SiC grains resulted in very low mass fractions of such 
grains, typically less than 1 per cent. This means that the 
conditions in the SMC (ISRF intensity, duration of the 
exposure by the UV radiation) are probably favourable 
for converting icy mantles into organic refractory, but not 
for the further processing of organic refractory into amor- 
phous carbon. In contrast to the SMC, the presense of icy 
mantles may be allowed for th e dus t grains in the Galactic 
diffuse medium (Zubko et al 



1998) 



Following Zubko et al. (1998), the models based on the 
silicate core-organic refractory mantle (composite) grains 
are refered to as G (M) models. The GM model is a combi- 
nation of G and M models and contains both core-mantle 
and composite grains. As shown in Fig. [I] all the above 
models fit the extinction curve quite well. The size distri- 
butions of both core-mantle and composite grains are quite 
wide and cover both small and larger grains with the pref- 
erence to the grains of sizes 10-100 nm. Silicon nanopar- 
ticles have a diameter of 3.0 nm by definition. All the 
models consume the maximum amounts of carbon, oxygen 
and silicon allowed for dust and slightly less for magnesium 
and iron. All the carbon consumed is contained in organic 
refractory Approximately equal amounts of silicon are 
locked up in the silicon nanoparticles (4-5 p.p.m.) and 
in other components (5-6 p.p.m.). Silicate core-organic 
refractory mantle grains prevail by mass in all the cases. 
Note that the dust composition in our models is drastically 
different from that in previous SMC models, which were 
based on modification of the standard MRN model (Pei 



ZUBKO 



1992 ; Rodrigues et al. 1997 ). The grain size distributions 
in our models are not a simple power law, as assumed in 
the previous models, but are rather optimized to repro- 
duce the extinction curve and to simultaneously obey the 
abundance constraints. 

Since we do not know presently the wavelength- 
dependent intensity of the ISRF in the SMC and, in addi- 
tion, the existing extinction curves for the SMC have a UV 
boundary at around 0.13 /itm, we are unable to estimate 
the fraction of the UV photons absorbed b y the silicon 
nanoparticles. In the Galaxy, Gordon et al. (1998) found 
this fraction to be 0.10±0.03. However, as was found by 



Zubko et al. ( 1999 ), the mass fraction of silicon nanoparti- 
cles may serve as a good indicator in this case. The values 
of this quantity derived for the models presented above: 
0.07, 0.071, and 0.085 suggest that the SMC is similar to 
the Galaxy, and therefore should also be a source of signifi- 
cant ERE. Moreover, the proximity of the silicon nanopar- 
ticle mass fractions obtained by modeling rather different 
extinction laws: Galactic and SMC, with different chemi- 
cal constraints and different dust models suggests that sili- 
con nanoparticles and ERE may be a universal phenomena 
in galaxies. 

It is evident from Fig. [I] and Table [l that each of the 
models presented is almost equally good (as indicated by 
p,) in fulfilling the requirements formulated above, with 
the G model being slightly more preferable. In order to 
discriminate between the models, we calculated the model 
scattering properties, albedo and asymmetry parameter, 
which are displayed in Fig. g. We compare our results with 
the obs ervati onal data for the Galaxy taken from Gordon 
et al. (1997) to illustrate the significantly different pre- 



dictions for the SMC. The model albedos are quite close 
to one other for all wavelengths, whereas the model asym- 
metry parameter shows large differences, especially in the 
UV. In general, the model albedos are lower than respec- 
tive Galactic ones, except for the near IR, where we may 
see an opposite effect. Only the asymmetry parameter of 
M model is close to the re spect ive Galactic values. Note 
that the albedos of Pei's ( |1992| ) model of the SMC dust 
significantly exceed the albedos, calculated in our models, 
and also the expected Galactic values (Fig. ||). Another 
possible mean to choose the most appropriate model may 
be the polarization data. We hope to include these data 
in a self-consistent analysis in forthcoming papers. 

In summary, we report here for the first time more re- 
fined models of the SMC bar dust which are in good agree- 
ment with the observed extinction, elemental abundances, 
and the strength of the ISRF. The models were calculated 
by using the regularization approach. The major grain 
constituents were found to be silicates, organic refractory 
and nanosized silicon. This conclusion is subject to some 
uncertainty due to the uncertain element depletion pat- 
terns in the SMC. We predicted the scattering proper- 
ties of our models, which are significantly different from 
the Galactic values. More observational constraints, e.g. 
the polarization data, are to be included into analysis to 
choose the most appropriate dust model. 

I thank Karl Gordon and Geoff Clayton for providing 
me the extinction curve for AzV 398 in the electronic ta- 
ble. During my work, I benefited from many stimulating 
discussions with Ari Laor. This research was supported 
by a grant from the Israel Science Foundation. 
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Fig. 1. — Size distributions of dust grains (left panels) and respective extinction curves (right panels), fitting the extinction curve toward 
AzV 398. Shown are G, M and GM dust models. The uncertainties of the size distributions were calculated using the method of statistical 
modeling by Zubko (in preparation). 
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Fig. 2. — Scattering properties: albedo and asymmetry parameter, corresponding to various dust models fitting the extinction curve toward 
AzV 398. The diamonds depict the observational data for the Galaxy taken from Gordon et al. (1997) and the circles the albedo of Pei's (1992) 
model of the SMC dust. 
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Table 1 

The main parameters of the models fitting the extinction toward AzV 398. 



Model 


Components 


C 


Si 


O 


Mg 


Fe 


/mass 


M 


A* 


SMC 


gas + dust 


46 


10 


107 


9 


7 










dust 


20 


10 


40 


9 


7 








G 




20 


10 


40 


5 


5 


100.0 


2.62-27 


0.211 




Silicate(0.4)-ORR(0.6) 


20 


5 


40 


5 


5 


91.5 








Si(0.99)-SiO 2 (0.01) 





5 











8.5 






M 




20 


10 


40 


6 


6 


100.0 


2.81-27 


0.269 




Composites (45%-porous) 


20 


6 


40 


6 


6 


93.0 








Si(0.99)-SiO 2 (0.01) 





4 











7.0 






GM 




20 


10 


40 


6 


6 


100.0 


2.75-27 


0.233 




Silicate(0.5)-ORR(0.5) 


13 


5 


31 


5 


5 


70.4 








Composites (45%-porous) 


7 


1 


9 


1 


1 


22.5 








Si(0.99)-SiO 2 (0.01) 





4 
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a First two lines contain the SMC elemental abundances in gas+dust and dust. Next lines present the data for the 
models reported: elemental abundances in ppm (atoms per 10 6 H atoms), mass fractions of each dust component (/ ma ss)> 
the total mass of dust matter in g per H atom (M) and the values of p, characterizing a quality of a model (see Zubko et 
al. 1996 for more details). The volume fractions of grain constituents are noted in brackets. The 45%-porous composite 
grains mostly contain silicate (0.28 in model M and 0.23 in model GM) and organic refractory residue (ORR)(0.24 and 
0.23) with less amount of amorphous carbon of various forms (0.09 and 0.03). 



